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Abstract. We introduce a family of 3D combinatorial objects that we define as minimal 3D polyominoes 
inscribed in a rectanglar prism. These objects are connected sets of unitary cubic cells inscribed in a 
given rectangular prism and of minimal volume under this condition. They extend the concept of 2D 
polyominoes inscribed in a rectangle defined in a previous work. Using their geometric structure and 
elementary combinatorial arguments, we construct generating functions of minimal 3D polyominoes in the 
form of rational functions. Wc also obtain a number of exact formulas and recurrences for sub- families of 
those polyominoes. 



1. Introduction 

Since the rise of modern combinatorics in the early 1960's, most combinatorial objects are visualized 
and investigated with pencil and paper and therefore, arc 2-dimensional. A number of extensions from 2D 
combinatorial objects to 3D objects were introduced: Ferrers diagrams were extended to plane partitions, 
permutations were extended to maps on a surface and to braids, 2D fractals were extended to 31? fractals 
and a short list of exact results for the enumeration of 3D objects have been produced so far (see [T],[S|). 
Behind these efforts lay a fundamental question: Is 3D combinatorics similar to 2D combinatorics in the 
sense that it is a natural extension of notions and concepts already known in 2D or does it introduce new 
material and concepts unknown in 2D combinatorics ? This question was part of our motivation to begin a 
study of 3D polyominoes. 

A 2D-polyomino is a 4-connected set of unit square cells in the discrete plane. That is, the cells are 
connected by their edges. A polyomino is inscribed in a b x k rectangle when it is contained in this rectangle 
and touches each of its four sides. Inscribed 2D polyominoes with minimal area were introduced in a previous 
work (see [1]) where an elementary geometric characterization was given that permitted their enumeration 
and the construction of their generating functions. The geometry of an inscribed minimal 2D polyomino 
can be described in simple terms as a hook- stair-hook structure where a hook is formed with two mutually 
perpendicular rows of cells starting on an edge of the rectangle and meeting at their corner end (see fig. 



3(a) red cells and [I] for more details). A 2D stair is a path of connected cells beginning on one corner of 



a rectangle, say north-west, and moving along the corresponding diagonal in the east - south direction (see 



fig. 3(a) black cells and their circumscribed rectangle) to end in the opposite corner of the rectangle. 

In this paper we introduce 3D polyominoes inscribed in a rectangular bxkxh prism and we define them as 
collections of unit 6-connected cubic cells contained in the prism and touching each of its six faces. We give a 
geometric description of a complete collection of families of inscribed 3D polyominoes with minimal volume. 
This allows us to present generating functions, recurrences and exact formulas for families of minimal 3D 
inscribed polyominoes. 

3D polyominoes, sometimes called polycubes, are known in the litterature and in recreational mathematics 
in the context of packing problems (see [2]) and their enumeration according to their volume is given up to 
volume 16 in [5] as the result of a computer program. However combinatorial enumeration of inscribed 3D 
polyominoes does not seem to have been considered so far. 

We will introduce three disjoint families of minimal 3D inscribed polyominoes and show that their union 
forms the complete set of 3D inscribed polyominoes with minimal volume. These three families will be called 
respectively 3D diagonal poliominoes, 2D x 2D polyominoes and skew cross polyominoes. 

We will use the orthogonal projection of inscribed 3D polyominoes on the upper face of the prism in view 
of the fact that an inscribed 3D polyomino is of minimal volume if and only if its orthogonal projection on 
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(a) Generic 2D minimal (b) 3D diagonal polyomino 

polyomino 

FIGURE 1. 2D and 3D inscribed minimal polyominocs 

each rectangular face of the circumscribed prism is a 2D polyomino of minimal area. This is easily proved by 
contradiction for if a 3D inscribed polyomino is not minimal, then one of its projections is not 2D minimal. 
Similarly, if one projection is not minimal, then the 3D polyomino cannot be minimal. 

The reader will observe that pictures illustrating 3D polyominoes frequently appear in our proofs. This 
visual support was used in our investigation and it helped us develop 3-D visualisation. Their absence would 
increase the difficulty of understanding our arguments. 

Notations. We will use capital letters for sets and generating functions and their corresponding lower case 
letters will be used for set cardinalities. For example P3D,min(P,k,h) will denote the set of 3D polyminoes 
inscribed in a b x k x h rectangular prism with minimal volume, P3D,min 

(b, k, h) will be their number and 

PzD,min{x, y, z) = J^b k h-P3£>,mm(&, k, h)x b y h z h will be their generating function. We will use the convention 
that the edge of length b of the cube is along the x axis and similarly the lengths k, h are along the y and z 
axis respectively. 

The degree of a 3D cell c in a polyomino, denoted deg(c), is the number of cells having a face in contact 
with c and the degree of a 2D cell c is the number of cells with an edge contact with c. All polyominocs 
considered in this paper are 2D or 3D, always inscribed in a rectangle or a rectangular prism and of minimal 
area or volume. Therefore we will often omit to specify these caractersitics of polyominocs. We will use 
trinomial coefficients in their standard notation ( a ^ + c c ) ■ We refer the reader to [[J]] for results and definitions 
on 2D polyominoes. 

The paper is organized as follow. In section^ we introduce diagonal 3D polyominoes and the subfamilies 
needed for their geometric description. We give generating functions, recurrences and exact formulas for 
a number of these subfamilies. In section [3j we define two families of non diagonal polyominoes: 2d x 2d 
polyominoes and skew cross polyominoes with other subfamilies necessary to their description. We give their 
generating functions and some exact formulas. In section [4j we prove the main result of the paper which 
states that these three families of polyominoes form a complete set of 3D minimal inscribed polyominoes. 
This result comes with a rational form for the generating function PiD,min{x, y, z) of minimal 3D inscribed 
polyominocs. 



2. Diagonal polyominoes 

In similarity with 2D stairs, we define a 3D stair as an inscribed polyomino of minimal volume forming 
a path starting in a given corner of the prism, say the north-west-back corner, and moving with unit steps 



in the south, east or forward direction until it reaches the opposite 3D diagonal corner as in figure 2(d) In 
what follows, we will use 3D stairs as components of polyominocs. 

Recall that a 2D corncr-polyomino is a 2D minimal polyomino inscribed in a rectangle with a cell in a 
given corner of the rectangle. The number P c (b, k) of 2D corner-polyominoes inscribed in a b x k rectangle 
satisfies the following recurrence and exact formula: 
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fb + k- 2\ 

(2) =2 , , -1 



(1) P c (b, k) = l + P c (b, k-l) + P c (b - 1, k) 

6-1 

with initial conditions P c {b, 1) = P c (l,k) = 1. Its generating function has the rational form 

(3) P c (x, y)=J2 P ^ k >"y k = n \ 7i ~W\ T 

b j^ 1 (1-x-y) (l-x)(l-y) 

Recall also (see that the total number of polyominoes of minimal area inscribed in a rectangle P2D, min (b, k) 
of size b x k is given by the formula 

P2D, min (6, fc) = 8 ^ | * ~ ^ + 2(6 + fc) - 3bk - 8 

We will first define and investiguate 3D corner-polyominoes. A 3D corner-polyomino is a minimal poly- 
omino inscribed in a prism with one cell in a given corner of the prism, say the north-west-back corner. Let 
P c (b, k, h) be the set of corner-polyominoes inscribed in a 6 x k x h prism. 

Theorem 1. For all positive integers b,k,h, the number p c (b,k,h) of 3D polyominoes inscribed in a prism 
of size b x k x h with minimal volume and one cell in a given corner of the prism satisfies the following 
recurrence : 



[ ^{bUt^Th-i) - 1 V b=1 or k=1 or h=1 

Pc(b, k,h) = ll + 2( 6 !Y) + 2{ b \ h _- 2 ) + 2( k+ k h _- 2 ) - 6 

[ + p c (b — 1, k, h) + p c (b, k — 1, h) + p c (b, k, h — 1) otherwise 

Proof. The first case is the 2D case. It provides the initial conditions for the 3D case and is obtained 
from equations ([TJ and ([2]). In the second case, observe that a corner cell has degree one, two or three. 
There is exactly one 3D corner-polyomino of degree three inscribed in a 6 x k x h prism and we call this 
polyomino a tripod. This explains the term 1 in the recurrence. When the corner cell c is of degree two, 
then c is the corner cell of a 2D corner-polyomino different from a 2D hook that is inscribed in a face of the 
prism and attached to a perpendicular row of cells along an edge of the prism. A row of cells connecting the 
polyomino to a face of the prism will often be considered and we will call these components pilars. Figure 



2(b) illustrates this situation: the corner cell of degree two is the red cell, the 2D corner-polyomino is made 
of the red and blue cells and the set of green cells forms a pilar. The next four terms in the recurrence are 
thus deduced from equation ((2|) 



Now if the corner c has degree one, as in figure 2(c) then the polyomino starts with a 3D stair giving the 



last three terms of the recurrence and the proof is complete. Observe that the separation according to the 
degree of the corner cell also gives the following equivalent formulation for the recurrence: 



D c (6, fc, h) = tripod + (2D-corner — 2D-hook) + degl 

= 1 + (P c (b, k, 1) + P c (b, 1, h) + P c (l, k, h) - 3) + 
(P c (b, k,h-l) + P c (b, k -l,h) + P c (b - 1, k, h) 



Generating function. To establish the generating function for the set of 3D corner-polyominoes, we will 
first give the generating functions Stair(x, y, z), Tripod(x, y, z), 2dhook(x, y, z) and Deg2{x, y, z) which are 
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(a) deg3: tripod (b) Degree two: 3-D (c) Degree one (d) 3D stair 

hook 



Figure 2. Corner polyominoes 



3D stairs, tripods, 2D hooks and 31? corner-polyominoes of degree two respectively: 
Tripod(x, y, z) 



J2 xl y 3zk 



2 2 2 

xyz 



(4) 



Stair{x, y, z) 



2Dhook(x, y, z) 

Deg2(x,y,z) 
2xz 



i,j,k>2 

E 

i,3,k>\ 



(l-x)(l-y)(l-z) 



i+j+k-3 
i - l,j - l,k- 1 



x l y 3 z k = xyz ^^(x + y + z)' 

n>0 
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(l-i) 

2xy 
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y) (l-x)(l-y) 



(1-z) 



(1-2/) 

The proof for the rational form of these generating functions is straightforward once we understand the 
geometric nature of the corresponding objects: there is one tripod per prism because, by definition, their 
corner cell is in a given corner of the prism. The number of stairs from one corner to its diagonal opposite 
corner in a prism of size b x k x h is equal to the trinomial coefficient (t_i'i'iii)' 2D-hooks appear on 
a slice parallel to one of the faces so we have three terms, one for each coordinate plane. The generating 
function for corner-polyominoes of degree two (equation (|4])) is directly obtained from its definition: a 2D 
corner of degree one perpendicular to a pilar. 

Now we are ready to use these building blocs. For instance a 3D corner of degree one always begins as a 
3D stair of length at least two connected to a 3D corner of any degree. The generating function Degl(x, y, z) 
of 3D corners of degree one is thus 



Degl(x, y, z) = (Stair (x, y, z) — xyz) (1 + Tripod + Deg2 + 2Dhook) 

Since we now have the generating functions for corner-polyominoes of degree one, two and three, we 
deduce the following result. 

Proposition 1. The generating function P c (x,y,z) for 3D corner-polyominoes is the following : 
P c (x,y,z)= Pc(b 1 k,h)x b y k z h 



(5) 



b,k,h>l 

Stair(x, y, z) 



^ ^ Tripod(x, y, z) + Deg2(x, y, z) + 2Dhook(x, y, z) 

xyz 



Proof. This is an immediate consequence of the fact that a 3D corner-polyomino is the connection of a 
3D stair with a 3D corner-polyomino of arbitrary degree. In equation ([5]), we decide that the corner cell 

4 



common to a 31? corner-polyomino and a 3D stair belongs to the stair so we divide the generating function 
of the former by xyz. 



Theorem 2. For all positive integers b,k,h, we have 



(6) 



p c (b,k,h) = 4 
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h - 1 
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b + h-A-2i\ (b+k+h-A-i 



h-2-i 



b + h-3-2i 



+ 



h- 1 
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Proof. By induction on b + k + h. li h = \ then the prism is reduced to a rectangle in the xy plane and 
formula © gives p c 

(b,k,l) = 2( & |V) - 1 which 

agrees with equation @. The same argument is true for 
b = 1 and fc = 1. Suppose that formula ([3]) is true for a prism of size b x k x h with b,k,h > 2. We have 



p c (6,fc,/i+l) = 1 + 2 



"6 + k - 2\ /& + h - 1\ fk + h-1 

+ p c (b- l,k,h + l)+p c (b,k- l,h+l) +p c (b,k,h) 
by theorem [1] and by induction hypothesis we obtain 
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p c (b,k,h + l) = 1 + 2 
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(l + (-lf +1 ) 
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(l + (-lf +1 ) 



= 4 



6 + h - 1 

6-1 



6 + h - 1\ /6 + k + h - 2 
b + h-l 

b + k-2 
fc-1 



'fc + /i - 2 
/i - 1 

h-l 



b + h- A-2i\ fb + k + h- A-i 



h-2-i 



b + h-3-2i 
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6 + /i-3-2A/6 + fc + /i-3-i 



h-1- 



(l + (-lf +1 ) 



6 + - 2 - 2z 



which is what we wanted to prove. 

It is now possible to construct formulas for the set of polyominoes along one given diagonal of the prism. 
We define diagonal polyominoes as inscribed polyominoes of minimal volume formed with three pieces: two 
hooks on each corner of a diagonal of the prism connected by a stair in contact with their corner cell (see 
figure 1(b) ). By a hook we mean either a 3D corner-polyomino with corner of degree two or three or a 2D 
hook. From this definition, we deduce the rational form of the generating function of diagonal polyominoes. 

Proposition 2. The generating function lDiag(x,y, z) of diagonal polyominoes along one given diagonal 
of a prism is the following 
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(7) 



lDiag(x, y, z) = Stair(x, y, z) 



1 



Tripod(x, y, z) + Deg2(x, y, z) + 2Dhook(x, y, z) 
xyz 
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(a) 3D cross 



(b) A polyomino on two diagonals 



Figure 3. Diagonal polyominoes on more than one diagonal 

Proof. This is a direct consequence of the definition of diagonal polyominoes, tripods, stairs, corner- 
polyominocs and 2D hooks. Notice that the number 1 inside the brackets of equation (J7J stands for the fact 
that 3D hooks could be absent and we divide by xyz the next term because we arbitrarily decide that the 
cell common to a hook and a stair belongs to the stair so that we remove it from the hook with this division. 

In the next step, we want to count the total number of diagonal polyominoes in a prism. There are four 3D 
diagonals in a prism. If a polyomino belongs to exactly two diagonals, then the two diagonals always define a 
plane perpendicular to two parallel faces of the prism. The orthogonal projection of the polyomino on these 
faces must be a 2D minimal polyomino and therefore this projection has the generic form hook- stair-hook of 
a 2D minimal polyomino. The projection of the two 3-D diagonals on any other face are the two diagonal of 
these rectangles. Since the only 2D polyomino that belongs to two diagonals of a rectangle is a 2D cross, i.e. 
two perpendicular rows of cells, the projection of the polyomino on the other faces is always a 2-D cross. This 
has consequences on the form of any 3D polyomino along two diagonals which must be made of a full pilar, 
i.e. a pilar connecting two opposite faces, connected to a perpendicular 2D generic polyomino inscribed in a 
full 2D slice of the prism (see the blue part in figure [3(b)] Moreover the full pilar must meet the orthogonal 
2D polyomino on its stair part. 

Now if a diagonal polyomino belongs to three diagonals, then its projection on each face of the prism is 
a 2-D cross. The only polyomino whose projection on all faces is a 2D cross must be a 3D cross which also 



belongs to four diagonals (see figure 3(a) ). 



Polyominoes along two diagonals. The generic form of polyominoes on two diagonals can be described 
as two 2D corner-polyominoes sharing their corner cell which belong to a full pilar perpendicular to the 
corner-polyominoes. Since we already know the generating function for 2D corner-polyominoes, it is easy to 
deduce the generating function for diagonal polyominoes belonging to two and three diagonals. 

Proposition 3. The number 2diag z (b,k,h) of 3D diagonal polyominoes belonging to the two diagonals 
perpendicular to the xy face of a prism such that the projection of these two diagonals has a vertex in the 
upper left corner of the face of size b x k has the following generating function 

2diag z (x,y, z) = 2diag z (b,k,h)x b y k z h 



b,k,h>l 

2 

x 



1 / 2xy xy 



xy\(l-x-y) (l-x)(l-y)J (1 - zf 

Proof. This is immediate from equation Q and the fact that these polyominoes have the geometric 
structure 2D comer x (2D corner — corner cell) x pilar. 

3D crosses. Next we need the generating function 3Dcross(x, y, z) of 3-D crosses which are the 3-D minimal 



polyominoes made only of pilars, at least three, meeting on one common cell c (see figure 3(a)). Observe 
that for a prism of size b x k x h with b,k,h > 2, there are bkh cross polyominoes inscribed in that prism 
and only one if any two of these three parameters equals one. We will only consider crosses in a box of size 
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at least 2x2x2. We thus have: 

, , , h h h x 2 (2 ~ x) y 2 (2 - y) z 1 (2 - z) 
3Dcross(x,y,z) = > bkhx b y k z h = — i %f 1 — ^ 

Proposition 4. The generating junction Diag(x,y, z) of the total number of diagonal polyominoes is the 
following 

Diag{x,y,z) = diag(b,k,h)x b y k z h 

b,k,h>2 

(8) =4x lDiag(x, y, z) — 2 x {2diag z (x, y, z) + 2diag y (x, y, z) + 2diag x (x, y, z)) + 3 x 3Dcross(x, y, z) 

Proof. In order to count all 3D diagonal polyominoes, we use some inclusion-exclusion principle. Here 
are the steps: 1- Count polyominoes along one diagonal and multiply by four. 2- The polyominoes that 
belong to two diagonals or more were counted twice or more so for each pair of 3D diagonals, remove the 
polyominoes belonging to those two diagonals. 3- The polyominoes belonging to three diagonals, and thus 
to four, were counted four times in the first step, removed six times in the second step and so must be added 
three times to be counted once. 

Notice that this inclusion-exclusion argument is not valid for degenerate prisms that have one side of 
length one and for their corresponding terms in the generating function ([5]). 

In table [T] below, we present the numbers diag(n,n^n) of 3D diagonal polyominoes as well as cardinality 
of other sets of polyominoes that are inscribed in a cubic prism of size n x n x n that are described in this 
paper. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


diag(n, n, n) 


1 


32 


2271 


79936 


2103269 


49998072 


1163531779 


27263453288 


P2Dx2D(n,n,n) 








66 


2256 


34092 


352992 


2994750 


22756896 


SC a (n,n,n) 








48 


3456 


85008 


1321344 


16174416 


172476672 


SCt(n, n, n) 








16 


1408 


33776 


505472 


5998512 


62474496 


P3D,min{n,n,n) 


1 


32 


2401 


87056 


2256145 


52177880 


1188699457 


27521161352 



Table 1 . Minimal polyominoes inscribed in a n x n x n prism 



3. NON DIAGONAL POLYOMINOES 

Does there exists minimal 31? polyominoes that are not diagonals ? The answer is yes and figure U shows 



a sample of these objects. For instance the polyomino in figure 4(a) is not diagonal because it is impossible 
to identify a corner-polyomino in a subprism having three faces in common with the circumscribed prism. 
This polyomino can be seen as the juxtaposition of two perpendicular 2D polyominoes each with contact 
cell that is not a corner cell. This is our definition for the family of non diagonal minimal polyominoes that 
we call 2D x 2D polyominoes. 

3.1. 2D x 2D polyominoes. In what follows, we establish the generating function for 2D x 2D polyominoes. 
For that purpose, we split these polyominoes in three parts, each part corresponding to one color in figure 



4(a) The central part, made of green cells with red corners, will be called a skew hook. It consists of three 
mutually orthogonal segments of cells. The two end segments touch a face of the prism and so are pilars 
with at least one cell. They touch the middle segment on its end cells. These two end cells are the contact 



cells of the two other parts (one in blue and one in yellow in figure 4(a)). If we discard the two pilars, 
each end cell of the middle segment can be seen as the corner cell of a 2D corner-polyominoe. The two 2D 
corner-polyominoes with their associated pilars are perpendicular and each one goes from one face to its 
opposite face. Notice that the smallest prism that contains a 2D x 2D polyominoe has size 2x3x3 and in 
that case, the polyominoes are made of two perpendicular full pilars that are the discrete version of euclidian 
skew lines. 

We begin with the generating function of skew hooks. This is quite elementary when we consider that each 
pilar contains at least one green cell and the central segment contains two red end cells but not necessarily 
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(a) case 1 



(b) case 2 



(c) case 3 



Figure 4. Non diagonal Polyominocs 



green cells. In order to fix ideas, we agree that the yellow 2D polyomino is in the yz plane with z length at 
least two if we count the red corner cell. The blue 2D polyomino is in the xy plane. If we decide that we 
do not count the contribution in x and z of the central segment and the contribution in y of the red corner 
cells, we have the following generating function SH(x, y, z) for skew hooks : 

SH(x,y,z) 



(1-x) (l-j,) (1-z) 

The yellow 2D corner polyominoes in the yz plane of z height at least 2 are obtained from equation ([3]): 

/ \ / 2 1 1 
2-t> c , z >2 [y, z)=yz[ - — r - 

Then the blue 2D corner polyominoes in the xy plane of x length at least 2 : 

2 1 1 



2D c , x > 2 (x, y) = xy 



(l-x-y) (l-x)(l-y) (1-y) 



In order to assemble these three components, observe that if we fix the vertical pilar and the yellow 2D 
corner, then the horizontal green pilar may take two directions that determines the direction of the blue 2D 
polyomino which is equivalent to multiply by two the number of blue 2D corner-polyominoes and remove the 
2D crosses which would be counted twice otherwise. We do the same for the yellow 2D corner-polyominoes. 
Finally, observe that the yellow polyomino could be on the left rather than on the right of the central part 
which multiplies by two again the number of polyominoes and we obtain the following generating function 
for 2D x 2D polyominocs with orthogonal planes xy and yz. 

(9) P xyxyz (x, y, z) = 2 (2 ■ 2D CtX > 2 (x, y) ■ 5* ■ (2 ■ 2D c<z > 2 {y, z) (1 _ ^ _ - } 

Finally, observing that two pairs of orthogonal planes determine two disjoint sets of 2D x 2D polyominocs, 
we obtain the generating function P2Dx2d(x, y, z) for the total number of non diagonal 2D x 2D polyominoes 
by adding the three generating functions corresponding to each pair of orthogonal planes: 

(10) P2Dx2D(x,y,z) 

— Pxyxyz + Pxyxxz + Pxzxyz- 

3.2. Skew cross polyominoes. We define our second family of non diagonal polyominoes as follow. A 
skew cross polyomino starts with a central cell c of degree three which is the corner cell of three 2D corner- 
polyominoes mutually perpendicular. We partition this family in two types. Type a) The cell c has two 
parallel contact faces. Type b) The three contact faces of the central cell c are incident to a vertex of c. 
These two families are illustrated in figure [5] 



(a) (b) 

Figure 5. Skew cross polyominoes 



Type a). We start by establishing the generating function for each of the three 2D corner-polyominoes 
needed to obtain a skew cross polyomino of type a). To fix the ideas, suppose that the three contact faces 
of the cell c have already been chosen and that the 2D corner-polyomino red and green is in the yz plane, 



the yellow part is in the xz plane and the blue part is in the xy plane as illustrated in figure 5(a) We 
have the choice between the red central cell c and the cell in contact with it as the corner cell of the 2D 
corner-polyomino. We choose the cell in contact with c. For the 2D corner-polyomino in the yz plane, the 
z length must be at least 2 and the generating function is 

(11) P c ,z>2(y,z) = yz 
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(12) P c ,x>2{x, y) = xy 

(13) P c ,z>2(x, z) = xz 

The product of the three series (fTTj) . (fT2 j) . (fT3|) gives the generating function of skew cross polyominoes of 
type a with preselected faces of the central cell provided we adjust with the fact that the z length of the 
cell c was counted twice and its y lenght was not counted. Now once the faces of c are chosen, there is some 
freedom for the direction of the 2D corner-polyominoes. Indeed, if we choose first one of the two directions 
of the corner-polyomino coming from the face between opposite faces, then we still have to choose between 
two directions for another 2D corner-polyomino. For two faces in the xz plane, we have two choices for a 
face yz. Thus the generating function SC a i{x, y, z) of skew crosses of type a when two faces in plane xz and 
one face in the plane yz are chosen is : 

Ay 

SC al (x 7 y,z) = — P c ,z>2(y,z)P c . x > 2 {x,y)Pc,z>2{x,z) 
^ 4x 3 y 3 z 3 (l + x- z){l +y- z)(l + y - x) 



(1 - x) 2 (l - y) 2 (l - z)2(l - y - z){\ - y - x){l - x - z) 



Knowing that there arc 12 triplets of faces of type a on a cell c, we sum six generating functions similar to 
equation (| 14[) and obtain the generating function SC a (x, y, z) for skew crosses of type a which simplifies to : 

-16x 3 y 3 z 3 ((l -x + y)(l - x + z) + (1 - y + x)(l - y + z) + (1 - z + x){l -z + y)) 
qIX ' 2/ ' Z) (1 - x) 2 {l - y) 2 {\ - z) 2 (l -y- z){\ -x- y)(l - x - z) 

Type b. To establish the generating function of skew crosses of type b, we choose three faces of the cell c 
incident to one vertex of c. We choose each cell in contact with a face of c to be the corner cell of a 2D 
corner-polyomino. There are two possibilities once the three corner cells are chosen. Here is the generating 
function for a given set of three faces corresponding to one vertex of c : 

Pc,z>2(x, z) x P c ,x>2(x,y) x P c ,y>2(y,z) + P c , y > 2 (x,y) x P c , z >2(y,z) X P c ,x>2(x, z) 
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There are 8 sets of three faces of c incident to one vertex of c and for each of these sets, we obtain the same 
generating function which means that the generating function for skew crosses of type b is the following: 

SC b (x, y, z) = 8(P C:Z > 2 (x, z) x P c , x >2{x, y) x P c , y > 2 (?/, z) 

+ Pc,y>2{x, y) X P c , z > 2 (l/, z) X P c ,x>2{x, z)) 

16x 3 y 3 z 3 ((l -x + y)(l - x + z) + (1 - y + x)(l - y + z) + (1 - z + x)(l - z + y) - 4) 
~ (1 - .t) 2 (1 - y) 2 (l - z) 2 (l -y- z){\ -x- y)(l - x - z) ' 

To obtain the generating function for all skew crosses SC(x,y, z), we simply add the generating functions 
for types a and b and observe that they have the same denominator and that terms in the numerator cancel 
except for a constant so that we obtain 

6Ax 3 y 3 z 3 



(15) SC(x,y,z) 



(1-x- y)(l -x- z)(l -y- z)(l - x) 2 (l - y) 2 (l - z) 2 ' 
4. Main result 



So far we have established the generating function for three disjoint classes of 3D polyominoes. We claim 
that the union of these three classes forms the whole set of 3D inscribed polyominoes with minimal volume. 

Theorem 3. The total number P3D,min{b 1 k,h) of polyominoes inscribed in a b x k x h rectangular prism 
and minimal volume b + k + h — 2 is the sum of diagonal polyominoes and non diagonal polyominoes of type 
2D x 2D and skew crosses: 

P3D,min(b, k, h) = diag(b, k, h) +p2Dx2D{b, k, h) + sc(b, k, h). 

Proof. In order to prove this result, we introduce a second classification of 3D polyominoes and we show 
that every set of polyominoes forming this classification belongs to one of our three families of polyominoes. 

Consider the orthogonal projection II(P) of an inscribed 3D polyomino P on the upper face of the prism. 
Observe that II(P) is a 2D inscribed polyomino of minimal area and therefore possesses the geometric 
structure hook- stair-hook of minimal 2D polyominoes. Two cells of the 3D polyomino will play a special 
role in our classification. We call them contact cells and define them as follow. For every polyomino 
P G P3D,min{b-,k,h) there is a unique 3D stair connecting the lower and upper faces of the prism which 
forms a non decreasing path from floor to ceiling. The two contact cells c\ , c 2 arc respectively, the last 
cell touching the floor and the first cell touching the ceiling in this path. We will use the positions of the 
projections II(ci), II(c2) in our classification. If, without loss of generality, we fix a 2D diagonal in the upper 
face to give a direction to the hook- stair-hook structure, there are ten positions of the pair II(ci), II(c2) with 
respect to the upper hook, each pair giving a class in this classification of P3D,min(b, k, h). The ten positions 
can be seen in figure U] where II(ci) and n(c2) are in black. Observe that these ten cases do not form a 
complete partition of the set P S P3D,min(b, k, h) but our goal is to provide a complete set of representatives 
up to symmetry so that every other case is similar to one of the cases considered. 

Another geometric observation used throughout this proof is the fact that minimality forces the cells that 
are not part of the 3D stair from c% to c 2 to form connected components that are all horizontal and attached 
to the 3D stair. Therefore the only cells that do not appear in the projection II(P) belong to the 3D stair 
between the two contact cells. We will show that the polyominoes belonging to each of the 10 cases also 
belong to one of the three families of polyominoes, namely diagonal, 2D x 2D and skew crosses. 

Cases 1 and 2. The polyominoes of cases 1 and 2 are 2D x 2D polyominoes and are illustrated in figure 
4(a) and |4(b)] Indeed the fact that the contact cells are projected on the same arm of a 2D hook implies 



that the 3D stair is in a vertical plane. The other cells of the polyomino form two connected components. 
One component is a pilar that ends the hook arm containing II(ci),II(c 2 ) and the other contains the other 
hook arm which is a pilar and a 2D corner. 



Case 3. In case 3, there is a contact cell projected on each arm of a 2D hook (see figure 4(c) |. The part of 
the polyomino that is not projected on this hook is a 2D horizontal corncr-polyomino that is neither on the 
lower face nor on the upper face of the prism. The two parts projected on each arm of the hook form also 
a 2D corner-polyomino so that we have a total of three 2D corner-polyominocs mutually perpendicular all 
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Figure 6. Classifcation of the projection of 3D polyominoes on the upper face of the prism 







(a) 



(b) (c) 

Figure 7. Polyominos of cases 4, 5 and 6 



(d) 



connected to a central cell projected on the corner of the hook thus forming a skew cross. 

Case 4. In case 4, we have a full vertical pilar from floor to ceiling projected on the corner of a 2D hook. 
Thus each part of the polyomino projected on an arm of the hook must be a horizontal pilar and so the re- 
maining part of the 3D polyomino must be a 2D horizontal corner-polyomino. We thus have three horizontal 
connected components at various heigths with two possible scenarios. Either the 2D corner-polyomino is 
between the two pilars (figure 7(a) I or it is not (figure [7(b) | . In the first case, we have a skew cross. In the 
second case we have a 3D diagonal polyomino. 



Case 5. In case 5 (see figure 7(c) I, there is a 31? stair in a vertical plane which is not a pilar and is pro- 
jected on the arm of a 2D hook, two horizontal pilars whose projection completes the arms of the hook. 
The remaining component must be a 2D horizontal corner-polyomino. There are again two subcases. Ei- 
ther the 2D corner-polyomino is at a height between the height of the two horizontal pilars either it is not. 



In the first case, the polyomino is a skew cross (figure 7(c) ). In the other case, the polyomino is a 3D diagonal 



Case 6. Here we have a 3D stair between the two contact cells whose projection joins a cell inside a hook 
with a cell outside the hook. Again, since the other connected components are horizontal 2D polyominoes, 
there are two horizontal pilars whose projection complete the arms of the hook. These constraints are suffi- 
cient to imply that the polyomino is a 3D diagonal. 



li 



(a) (b) (c) (d) 

Figure 8. Polyominos of cases 7, 8, 9, 10 



Case 7. Here the 3D stair connecting the two contact cells is projected on the 2D stair so that the remaining 
cells form 4 horizontal pilars, or 3 pilars when there is no second hook on IT(P) (figure 8(a)). No matter 
what the relative heigth of these pilars is, they always form a 31? diagonal polyomino. 



Cases 8 and 9. In these cases, the polyomino can always be decomposed in three parts: one 3D stair 
connecting the contact cells and two 2D horizontal corner-polyominoes with corner in contact with the 31? 
stair (figures [8(b)1 and [8(c)| )). This always give a 3D diagonal polyomino. 



Case 10. Here the contact cells are projected on different 2D hooks. These polyominoes always contain a 
3D stair between the two cells projected on the corners of the 2D hooks. Each of these two cells is the corner 
cell of a 3D hook. Thus the three parts form a 3D diagonal (see figure [8(d)] ) polyomino. 



5. Exact formulas 

In theorem [3] and other results of this paper, we chose to break generating functions in several parts 
because their full explicit formulation is too long. We did not provide exact formulas corresponding to all 
generating functions for similar reasons: the exact expressions are not always reducible. For example, here is 
an exact expression for the number sc(b, k, h) of skew crosses inscribed in a b x k x h prism true for integers 
b > 3, fc > 3, h > 3 that we could not reduce. 



(16) sc(b,k,h) = 64 ]T E E 

i=0 r=0 j=Q 



b+k—6 i b—Z—T / 7 \ / , \ / 7 

b \ I k \ I n 



3 + t + j J \3 + i + j — r ) \3 + i 



but if we turn our interest to the number of all minimum inscribed polyominoes of a given volume n, we 
obtain interesting exact formulas that lead to asymptotic information. In what follows, we will give exact 
formulas for each of the three families of 3D polyominoes presented in sections 2 and 3 to obtain one for the 
set P3D,min(n) of inscribed minimal polyominoes of volume n. 

Skew crosses. The recipe is the same for the three families : setting x = y = z in the generating function 
(|15|) . wc obtain the generating function 

«?(*)-£«<»)-»»- a _„/(!_,). 

and the exact formula for the number sc(n) of skew crosses of volume n: 

, ^ , s ,,4.9 / 9 \ fn b lln 3 9 844n 

(17) sc(n - 2) = 2 n+2 (n 2 - 27n + 194) - 8 — + + 12n 2 + —— + 96 
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2D x 2D polyominoes. We set again x = y = z in equations and (fTU| to obtain 

n i \ <• \ n+2 6x^(1 + 2a;) 2 

n>l y ' K ' 

„ , 2 , , / 3 „ 33 - 65 , 183 „ 544 , 147 \ 
(18) P 2Dx2D(n - 2) =3 • 2"+ 2 (n - 15) + i-n 6 - -n 5 + -n 4 - — n 3 + — n 2 + — n + 234 J 

Diagonal polyominoes. For diagonal plyominoes, we have to modify the generating function Diag(x, y, z) 
so that it becomes exact also for terms containing one of the variables x, y, z with power 1. This modification 
is done by removing from lDiag(x, y, z) the degenerate cases counting 2D polyominoes which cannot be part 
of the inclusion-exclusion calculus in equation [5] Then putting x = y = z we get 

, s . . ti _lo x 3 (36a; 8 + 129a; 6 - 207a; 5 + 234a; 4 - 126a; 3 + 49a; 2 - lOx + l) 

Diag(x) = 2, diaginp 

n>3 



)x n+2 = 



(1 -3x)(l -2x) 2 (l -xf 



1 91 / K<J 1 e 

(19) diag(n 2) =—T - 2"(45n - 411) - ( — it> - 



823 



6n 2 



22711 



4995 



8 24 60 16 

so that, adding equations (|17[) . (j 18[) . (|19|) .wc finally obtain an exact formula for P3D,min{n). 

Proposition 5. The generating function and exact formula for the numbers P3D,min{n) of 3D inscribed 
minimal polyominoes of volume n are 

t (n)x"+ 2 



P 



3D,min 



(n) 



x 3 {72x 10 + 36a; 9 + 510a; 8 - 1117a; 7 + 1276a; 6 - 1155a; 5 + 710a; 4 - 293a; 3 + 81a; 2 - 13a; + 1) 



(l-3a-)(l-2a;) 3 (l-a;) 7 



P3D, 



ll 2 ■ 3" +1 

16 

3n 6 



+ 2" +2 (4n 2 - 125n + 741) 
9?i 5 7n 4 

To" ~ ~ 



31727n 47739 



20 



16 



n 
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2 


3 


4 
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6 


7 
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10 


P3D,min(n) 


1 


3 


15 


83 


450 


2295 


10834 


47175 


190407 


719243 



Table 2. Numbers P3D,min{n) of minimal polyominoes of volume n 



An exact formula for P3D.mm(2, b, k). It is possible and relatively easy to derive one variable formulas for 
the numbers P3D,min(2, b, k) for different values of b and k with Maple and gfun but in the next proposition 
we develop combinatorially a two variables expression for P3D,mm(2, b, k). 



Proposition 6. The numbers P%d 

(20) P3D,min(2,b,k) = 



l (2,b^k) satisfy the expression 

b + k-2 
b-l 

16(6 + k - 2) - 12bk] (b + k 



4(b + k) 



(26 + 2fc-3)+4(6-2)(fc-2) 
1) 



Proof. We use the fact that a 3D minimal polyomino inscribed in a 2 x b x k prism is obtained by adding 
one cell to a 2D minimal polyomino inscribed in a b x k rectangle so that these 3-D polyominos are obtained 
from 2D polyominos rooted on one cell when we set the weight of the rooted cell x to be 2 deg ^ . We have 



(21) 



P 3Dmin(2,b,k)= £ ]T 2 de ^ 

peP2L>min(2,b,k) x£p 

Most cells in a minimal 2D polyomino have degree 2 except the leaves of degree 1 and one or two cells of 
degree 3 or 4. In order to transform equation [2T1 into equation [20l we have to partition the set P2D,min(b, k) 
with respect to the number of leaves and the number of cells of degree 3 and 4 . 
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Two leaves. 2D Polyominoes with two leaves are the 2( b ^~ 2 ) stairs polyominoes inscribed in a 6 x k 
rectangle. In such a polyomino, there are 6 + k — 3 cells of degree two and two cells of degree one. The 
number of 31? polyominoes obtained from these polyominoes is 




Three leaves. Polyominos with three leaves are the 2D corner-polyominos minus the hooks and the stairs 
and there are 4( b ~j^~ 2 ) — 2(6 + k) of them in a b x k rectangle. Each of these polyomino possesses b + k — 5 
cells of degree two, three cells of degree one and one cell of degree three. The number of ZD polyominoes 
obtained from these polyominoes is thus 

f b + k - 2^ 
6-1 



(23) 



- 2(6+ k) 



[{b + k - 5) x 2 2 + 3 x 2 1 + 2 3 ] 



Four leaves. Polyominos with four leaves are either crosses or non degenerate hook-stair-hook structures. 
There are (6 - 2)(fc - 2) crosses with four leaves and 2( b + k ^ 2 ) + 4(6 + k - 2) - 36fc - (b - 2)(k - 2) non 
degenerate hook-stair-hooks in a 6 x k rectangle. In a cross there are 6 + k — 6 cells of degree two, four cells 
of degree one and one cell of degree four. In a non degenerate hook-stair-hook, there are 6 + k — 7 cells of 
degree two, four cells of degre one and two cells of degree three. This gives the following number of 3-D 
polyominoes : 

(24) (6 - 2){k - 2) [(6 + k - 6) x 2 2 + 4 x 2 1 + 2 4 ] + 

b + h k ~ + 4(6 + k - 2) - 36fc - (6 - 2)(fc - 2) [(6 + k - 7) x 2 2 + 4 x 2 1 + 2 x 2 3 ] . 

The sum of expressions (|22|) . (|23 |) and (|24|) gives equation (|20]). 

We can easily extend the argument in the above proof to obtain exact expressions for P3D,min(3,b,k) 
from two-rooted 2D polyominoes and so on for P3D,mm(4, b, k), etc. For P3D jm m(3, 6, k) with 6 > 2, k > 2, 
we obtain 

P3D,min@, 6, k) =8(6 3 + fc 3 ) - 12(6 3 fc + 6fc 3 ) - 246 2 fc 2 - 46(6 2 + k 2 ) + 41(6 2 fc + 6fc 2 ) - 93fc6 + 58(6 + k) - 8 
'b + k - 2^ 



6-1 



(46 + 4fc-l)(26 + 2fc-3). 



Remarks. 

(1) In parallel with this work, one of the authors (H. Cloutier), wrote two programs to count minimal 
inscribed polyominoes. One program uses formulas obtained from the projection IT(P) of the poly- 
omino on the ceiling of the prism. The other program runs through all 3D polyominoes and keeps 
only the needed ones. We used the datas obtained from these programs to validate the results of 
this paper. 

(2) The argument in the proof of proposition |6] can be used to obtain exact formulas for 3D polyominoes 
inscribed in a prism with volume min + 1 when we use the formula giving 2D polyominoes of area 
min +1 in ([!]). 

(3) In the introduction of this paper, we asked a question about the relation between 2D and 3-D 
combinatorics. In this work, we had surprises moving from 2-D to 3-D but the answer to the question 
is not clear to us yet and we postpone our judgement. 

(4) The diagonal subseries -P 3 _D, mm (i) = J2 n P3D, mm (n, n, n)t n obtained from PsD,min(x, y, z) by setting 
all equals the exponents of x,y,z satisfies a functional equation of degree six in P3D,min{t) with 
coefficients that are polynoms in t. But no exact expression for PsD,min{'n>,n,n) could be found. 
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